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The trace distance between two states of an open quantum system quantifies their distinguisha- 
bility, and for a fixed environmental state can increase above its initial value only in the presence of 
initial system-environment correlations. We provide experimental evidence of such a behavior. In 
our all-optical apparatus we exploit spontaneous parametric down conversion as a source of polar- 
ization entangled states, and a spatial light modulator to introduce in a general fashion correlations 
between the polarization and the momentum degrees of freedom, which act as environment. 
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I. INTRODUCTION 

The dynamics of an open quantum system S inter- 
acting with an environment E is usually described by 
means of a completely positive trace preserving (CPT) 
map on the state space of the open system [T] . The very 
existence of such a map generally requires that the ini- 
tial correlations between the open system and the envi- 
ronment can be neglected, i.e. Pse(0) = p s {0) ® p E (0) 
where p s (0) = Tr B {p SE (0)} and p B (0) = Tr s {p SE (0)}. 
However, this assumption is not always physically jus- 
tified, especially outside the weak coupling regime [2J. 
Therefore, different approaches to the description of the 
reduced system dynamics in the presence of initial corre- 
lations have been developed in recent years [SHE]. 

An approach for the study of initial correlations that 
is based on the use of the trace distance and that does 
not rely on the determination of any reduced dynamical 
map has been introduced in [5]. In particular, one can 
find a clear signature of initial system-environment cor- 
relations as follows: if the environmental state is fixed, 
the trace distance between any two reduced states can in- 
crease over its initial value only in the presence of initial 
correlations. 

Recently the open system dynamics of two qubits has 
been experimentally investigated in all-optical settings, 
where the system is represented by the polarization de- 
grees of freedom, and the environment by the spectral 
[TU] or by the momentum [IT] degrees of freedom. In the 
present paper we provide an experimental proof of the 
feasibility and effectiveness of the abovementioned theo- 
retical scheme for the detection of correlations, observ- 
ing the effect of initial system-environment correlations 
in the subsequent open system dynamics by means of 
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the trace distance. In particular, we show an increase 
of the trace distance between two reduced states, shar- 
ing the same initial environmental state, over its initial 
value on both short and long time scales. Despite the fact 
that a full tomographic analysis can be performed, thus 
showing that the experimental setup can cope with the 
most general situation, the growth of the trace distance 
can here be detected simply by exploiting visibility data, 
thus showing that the theoretical analysis can really lead 
to efficient detection schemes. 

The paper is structured as follows: In Section [TT] we 
briefly present the general theoretical scheme. In Section 
|III| we describe the experimental apparatus in some de- 
tails, whereas in Section |IV| we describe the use of the 
spatial light modulator to introduce system-environment 
correlations, and analyze the evolution of the trace dis- 
tance. In Section [V] we provide the full tomographic re- 
consturction of the state under investigation. Section |VI| 
closes the paper with some concluding remarks. 



II. UPPER BOUND ON TRACE DISTANCE 
EVOLUTION 

The trace distance between two quantum states p 1 and 



p is defined as 
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with Xk eigenvalues of the traceless operator p 1 — p 2 , and 
its physical meaning lies in the fact that it provides a 
measure for the distinguishability between two quantum 
states [12] ■ It is a metric on the space of physical states, 
so that for any pair of states p 1 and p 2 it holds < 
D(p 1 ,p 2 ) < 1. Every CPT map A is a contraction for 
this metric: D(Ap 1 , Ap 2 ) < D(p 1 ,p 2 ), a property which 
will be crucial in the following analysis. 

The dynamics of an open quantum system can be char- 
acterized by investigating the dynamics of the trace dis- 
tance between a pair of reduced states p\{t) and p 2 (t), 
which evolve from two different initial total states p\ E (0) 
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and p 2 SE (0) . The change in the distinguishability between 
two reduced states can be interpreted as an informa- 
tion flow between the open system and the environment 
[13] • Indeed, since the reduced states p k s (t) are obtained 
from the corresponding initial total states p k E (0) k = 1, 2 
through the composition of a unitary operation and the 
partial trace, the contractivity under CPT maps implies 
that 

D{p\{t), p 2 s (t)) - D(pl(0),p 2 s (0)) < J 12 (0) , (2) 

where 

j 12 (o) ee d(/4(o),pL(o)) - i>g£(o),/4(o)) . (3) 

That is, the increase of the trace distance during the time 
evolution is bounded from above by the quantity Ii2(0), 
which represents the information initially residing outside 
the open system [9j. It is important to notice that the 
bound Ii2(0) can also qualitatively reproduce non-trivial 
features of the trace distance dynamics even if it is far 
from being reached [H]. If the initial total states are 
uncorrelated and with the same environmental state, i.e. 
PL(0) = Ps(0) ® P E (0) and p 2 SE (0) = p 2 s (0) ® p E (0), then 
^12(0) = 0. Thus, for identical environmental states, one 
can find an increase of the trace distance 

D(pl(t),p 2 s (t)) > D(pl(0),pl(0)) 

at a time t only if some correlations are present in at 
least one of the two initial total states. 

III. EXPERIMENTAL SETUP 

In our all-optical experimental setup the total system 
under investigation consists in a two-photon state pro- 
duced by spontaneous parametric downconversion. We 
look at the evolution of the two-qubit polarization en- 
tangled state, which represents the reduced system, and 
trace out the momentum degrees of freedom, which are 
not observed and represent the environment. We exploit 
a programmable spatial light modulator (SLM) to im- 
pose an arbitrary polarization- and position-dependent 
phase-shift to the total state. A linear phase is set both 
on signal and idler beams in order to purify the state 
[15j . whereas an additional, generic, phase function may 
be imposed to introduce initial correlations between the 
polarization and the momentum degrees of freedom in a 
very general way. A further linear phase is then used as 
a time evolution parameter for the two-qubit state. 

The experimental setup is shown in Fig. [TJ A linearly 
polarized CW, 405 nm, diode laser (Newport LQC405- 
40P) passes through two cylindrical lenses which com- 
pensate beam astigmatism, then a spatial filter (SF) se- 
lects a Gaussian spatial profile and a telescopic system 
prepares a collimated beam with beam radius of 550 /im. 
A couple of 1mm Beta Barium Borate (S) crystals, cut 
for type-I down conversion, with optical axis aligned in 




SF 



FIG. 1: (Color online) Diagram of the experimental setup. 



perpendicular planes, are used as a source of couples of 
polarization and momentum entangled photons |161 117) . 
The process preserves the total energy and the transverse 
momentum. The half wave plate (H) set on the pump 
path rotates the pump polarization in order to balance 
the probability amplitudes of generating a |VV) couple 
of photons in the first crystal or an |HH) couple in the 
second one. The couples are generated around a central 
angle of ±3° and we select A = lOmrad with two slits 
set on signal (2) and idler (1) paths. Two long-pass fil- 
ter (F) with cut-on wavelength of 780 nm set behind the 
couplers are used to reduce the background and to se- 
lect about 60 nm around the central wavelength 810 nm, 
while the two polarizers (P) are used to perform visibil- 
ity measurements as explained later on. The delay time 
between the probability amplitudes of generating a |VV) 
couple in the first crystal or a |HH) couple in the second 
crystal reduces the purity of the state. A nonlinear crys- 
tal (DC) with the proper length and angle is set on the 
pump path and precompensates this temporal delay [l"8l - 
|2"2"] . At first order a linear position dependent phase shift 
on both channels between |HH) and |VV) photons arises 
from the angle dependent optical path followed by |W) 
photons which must traverse the second crystal [19]. 



IV. SYSTEM-ENVIRONMENT 
CORRELATIONS AND TRACE DISTANCE 
EVOLUTION 

A. Theoretical description of the experiment 

In our scheme the SLM performs two basic tasks. First, 
it allows us to engineer the initial state by the introduc- 
tion of an arbitrary phase f(0). Aside from this, it pro- 
vides the effective system-environment interaction term 
sensitive to both the polarization and the momentum de- 
grees of freedom, through the introduction of a linear 
phase a9, where a is the time evolution parameter. The 
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total system-environment state is thus given by: 
1 



V2 



d9d9' g{9)g{9') 



x (\ho)\h6') + e liae+f(e)) \ve)\ve') 



(4) 



The factorized form g(9)g{9') is justified by the large 
spectral distribution [11]. Moreover, g(9) is a Gaussian- 
like shape function with FWHM of 6 mrad. Because of 
the phase /(#), the state in Eq. Q is correlated, i.e. 

p SE (a) = \ip3E(a)}(ip 3E (a)\ ^ p s (a) <8> p E (a) , 

and this is true als for the initial total state, i.e. for a = 0. 
Upon tracing out the momentum degrees of freedom, the 
polarization state is given by 

p s (a) = l(\HH)(HH\+e(a)\VV)(HH\ 



■e*(a)\HH)(VV\ + \VV)(VV\), (5) 



(a) = / de\g{9)\ 2 e^ e+ f^ 



where 



Since the angular distribution g{9) is symmetric and we 
use odd functions f(9), the quantity e(a) is real and it 
equals the interferometric visibility V{a) = Re[e(a)]. 

In order to characterize the effect of the initial system- 
environment correlations via the trace distance, we have 
to monitor the evolution of two different polarization 
states obtained from two different initial total states hav- 
ing the same environmental state. We compare an ini- 
tially uncorrelated state p\ E {a), corresponding to Eq. Q 
for f(9) = 0, with an initially correlated state p 2 E (a) for 
a non-trivial function f(9). In this way, the reduced sys- 
tem states p|(cv) k = 1,2 are both of the form given by 
Eq. ([5]), with different efe(a). Note that the product state 
pl E (0) differs from p|(0) £g> p 2 E (0) only for an overall phase 
term in the integration over 9, which has no observable 
consequences on the dynamics of the polarization degrees 
of freedom. The trace distance between the two reduced 
states under investigation is then given by 



D(p 1 s {alp 2 s {a)) 



e 2 (a) I 



(6) 



d9\g(9)f 
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Different choices for the initial phase f{9) result in dif- 
ferent dynamical behavior of the trace distance. We have 
exploited this fact to analyze in detail the effect of ini- 
tial system-environment correlations on the subsequent 
evolution of the open system. 



B. Experimental results 

Experimentally, we have measured the quantity e(a) 
for f(9) = and f(9) — sin(A6>), exploiting its equal- 
ity with the visibility, obtained in the standard way by 



counting the coincidences with polarizers set at 45° , 45° 
and at 45°, -45° (see [20] for further details). The func- 
tions of the variable 9 are discretized by the SLM, and 
thus become functions of the pixel number n. The resolu- 
tion is given by h/D, where h = 100 /im is the pixel width 
and D = 330 mm is the SLM distance from the source. 
In our experiment the SLM introduces the functions 

4> x (n) = -a opt (n-ni)+b (7) 
(f> 2 (n, a) = a opt (n - n 2 ) + a(n - n 2 ) + f(n - n 2 ) , 

on the two beams respectively, where a opt — 0.1 
rad/pixel is an optimal slope used to achieve the maximal 
purification of the polarization entangled state, and the 
constant b is used to offset the residual constant term. 
The integers ni and n 2 are the central pixel numbers on 
the idler and on the signal beams. The experimental evo- 
lution parameter is then a — ah/D and is expressed in 
rad/pixel. 
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FIG. 2: (Color online) Trace distance and visibility as a func- 
tion of the experimental evolution parameter a, the two quan- 
tities are related through Eq. |6|. Full circles describe the 
trace distance between Ps(a), i.e. fin — 712) = 0, and p|(a) 
with f(n — TI2) — sin(A(n — ^2)), A = —0.6 rad/pixel. Full 
squares describe the trace distance between Ps(a) and p|(a) 
with f(n — 112) = r(n — n-z), t — 0.1 rad/pixel. Lines are a 
guide for the eye. Empty circles refer to visibility with the 
choice f(n — 712) = 0, whereas empty squares refer to the case 
in which initial correlations are introduced through the phase 
function f(n — 712) = sin[A(n — "-2)]- For the visibility the 
uncertainties are within the symbols. 

The trace distance is the quantity which reveals the 
presence and the effects of initial correlations, and its 
behavior is reported in Fig. [2] together with the visibil- 
ity which provides the raw data from which the trace 
distance can be extracted in the present case. In the fig- 
ure full circles describe the trace distance, as a function 
of the evolution parameter a, between the reduced state 
Ps(a) evolved from the initial total product state, i.e. 
f(n — n2) — 0, and the reduced state p 2 (a) related to the 
initial correlated state with /(n — n 2 ) = sin(A(n — n 2 )). 
The trace distance, after an initial decrease and a first 
small oscillation, presents a revival up to a value which 
is more than three times the initial one. As expected, 
the reduced system can access information which is ini- 
tially outside it, related to its initial correlations with 
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the environment. The trace distance reaches its maxi- 
mum around a — 0.6 rad/pixel, toward the end of the 
monitored time interval. The maximum of the trace dis- 
tance quantifies the total amount of information which 
can be accessed by means of measurements performed 
on the reduced system only [14] . Note that it can be 
shifted to smaller values of the evolution parameter a by 
decreasing the absolute value of A. Thus, by introducing 
a sinusoidal phase modulation via the SLM, we have ob- 
tained a behavior of the trace distance which highlights 
the presence of initial correlations and their effects in the 
subsequent evolution, also for long times [25] , 

The simplest choice for the phase f(n — n-z) in the 
initially correlated state p 2 SE (a) is a second linear phase 
aside from that containing the evolution parameter a, i.e. 
f(n — ri2) = t(ji — ni). Indeed, this corresponds to shift 
the initially uncorrelated state p\ E {a) forward in time by 
r. Then, from the visibility measurement, we can directly 
obtain the evolution of the trace distance between p\(a) 
and p 2 (a) with f(n — 112) — r(n — 712)- This is repre- 
sented by full squares in Fig. [2j for r = 0.1 rad/pixel. In 
this case the growth of the distinguishability between the 
two reduced states starts from the very beginning of the 
dynamics. As expected, the trace distance increases over 
its initial value, reaching its maximum value at a = 0.1 
rad/pixel and decreasing afterwards. The subsequent os- 
cillations can be traced back to the finite pixel size. No- 
tice also that by using a linear term, we cannot obtain 
a revival of the trace distance (as in the previous case) 
over its initial value for high values of a. Since in this 
case p 2 (a) — pl(a + T), the full squares in Fig. [2] also de- 
scribe the evolution of the trace distance between a pair 
of reduced states occurring at two different points, sepa- 
rated by r, of the same dynamics starting from the ini- 
tial total product state given by Pg E (0). From this point 
of view, the increase over the initial value of the trace 
distance indicates that the single evolution under inves- 
tigation is not compatible with a description through 
a dynamical semigroup At, which could be introduced, 
e.g., on the basis of some phenomenological ansatz. In- 
deed, the semigroup property A t+T = A t A r , together 
with the trace distance contractivity under CPT maps, 
would imply D(pl(t) :P l(t)) = D(A t pl (0), A t pJ (r)) < 
D(pI(0),pI(t)) = D(pl(0),pl(0)). However, in general 
one cannot discriminate in this way whether the devi- 
ations from the semigroup dynamics are due to corre- 
lations in the initial total state or to other sources of 
non-Markovianity [24] . 



maximum is taken over all the projectors II or, equiv- 
alently, over all the positive operators II < 1. Upon 
considering the subspace spanned by {|HH),|VV}} and 
the corresponding a x Pauli matrix, the maximum is here 
obtained from the projectors on the eigenvectors of a x , 
which indeed give back half the difference between the 
visibilities. However, in more general situations one could 
need a full tomographic reconstruction of the reduced 
states. This would be the case in the presence of non- 
real coefficients efc(a) or when dealing with partially or 
fully unknown states. For this reason, we have also per- 
formed state reconstruction by polarization qubit tomog- 
raphy. By means of a quarter-wave plate, a half-wave 
plate and a polarizer, we measure a suitable set of inde- 
pendent two-qubit projectors [251 26] and then use the 
maximum-likelihood reconstruction of the two-qubit po- 
larization density matrix. In Fig. [3] (left) we show the 
tomographic reconstruction of the polarization state just 
after the purification and without any initial correlation, 
i.e. for f(n — 712) — and a = 0. The visibility is 

0. 914 ± 0.006 (not exactly one mostly because of the 
large spectrum detected). In Fig. [3] (right) we report 
the two-qubit tomography for the state characterizing 
the maximum revival of the visibility in the presence of 
initial correlations given by f(n — n-i) = sin[A(n — 712)], 

1. e. at a — 0.6 rad/pixel. The corresponding visibility is 
0.605 ± 0.007. 




FIG. 3: (Color online) Tomographic reconstruction of the two- 
qubit density matrix just after the purification (left), without 
any initial phase, i.e. for f(n — 712) = and a — 0. The 
visibility is 0.914 ± 0.006. Tomographic reconstruction for 
f(n — 712) = sin(A(n — 712)) at a = 0.6 (right), i.e. at the 
maximum of the visibility revival [compare with Fig. [2] . The 
corresponding visibility is 0.605 ± 0.007 



VI. CONCLUSIONS 



V. STATE RECONSTRUCTION 

In order to reconstruct the trace distance evolution, 
we only had to perform visibility measurements to ac- 
cess the off-diagonal values Ci(a). From a mathematical 
point of view, this corresponds to explicitly determine 
the projector operator defining the trace distance via the 
relation D(p 1 ,p 2 ) — maxn Tr {n (p 1 - p 2 ) }, where the 



We have reported an experimental observation of the 
effect of initial correlations between an open quantum 
system and its environment by means of the trace dis- 
tance. In particular, we have shown the increase of the 
distinguishability between two reduced states, sharing 
the same reduced environmental state, over its initial 
value on both short and long time scales. Our all-optical 
scheme is based on the use of a spatial light modulator, 
which allows us to introduce initial correlations in a very 



5 



general way. In particular, this setup allows to engineer 
different kinds of dynamical behavior of the trace dis- 
tance, so that one can, e.g., tune the position and the 
amplitude of the revival points of the distinguishability. 

Note added in proof. Recently, we became aware of 
[27], where initial correlations between the polarization 
and the spectral degrees of freedom of single photon 
states are experimentally witnessed by means of the trace 
distance. 
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